Anomalous oscillations of average transient lifetimes near crises 

Krzysztof KacperskiQ and Janusz A. HolystS 
Institute of Physics, Warsaw University of Technology, Koszykowa 75, PL-00-662 Warsaw, Poland 

and 

Max Planck Institute for Physics of Complex Systems, Nothmtzer Str. 38, 01187 Dersden, Germany 

(February 5, 2008) 



It is common that the average length of chaotic transients appearing as a consequence of crises 
in dynamical systems obeys a power low of scaling with the distance from the crisis point. It is, 
however, only a rough trend; in some cases considerable oscillations can be superimposed on it. In 
this letter we report anomalous oscillations due to the intertwined structure of basins of attraction. 
00 , We also present a simple geometrical model that gives an estimate of the period and amplitude of 

0^ ' these oscillations. The results obtained within the model coincide with those yielded by computer 

0^ | simulations of a kicked spin model and the Henon map. 
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I. INTRODUCTION 



Among a large diversity of phenomena investigated in the scope of deterministic chaos in nonlinear dynamical 
' systems one finds crises §}. These are sudden changes in the structure of a chaotic attractor due to collision with an 
C****" I unstable periodic orbit when a system parameter p crosses some critical value p c . 

In the dynamics after crisis (for p > p c ) characteristic transients appear; the system spends some time on the former 
(pre-crisis) attractor which is now a chaotic saddle. In the case of a boundary crisis the transient is followed by a 
definitive escape to some other attractor in the phase space while after an interior crisis transients are interrupted by 
(typically short) bursts to an extension of the pre-crisis attractor. After attractor merging crisis we have intermittent 
0^ , jumps between symmetric pre-crisis attractors. Such a behaviour is called crisis-induced chaos-chaos intermittency. 

For a large class of dynamical systems the time t that the system stays on the pre-crisis attractor has an exponential 
distribution Pr(t) = y exp(— y) with a mean value T fulfilling a power scaling low 

-a- 

£ T~(p-p c )-r. (1) 

For two-dimensional maps, the exponent 7 > can be expressed in terms of the eigenvalues of the periodic orbit 
q I involved in the crisis . The power low ([!]) and the formulas for the critical exponent 7 [||,[| have been confirmed 
r* ■ in many numerical and experimental studies (e.g. j|,[|). Throughout the recent years many other aspects of crises || 
. ! as wen as new types of crises has been investigated. 

However, it has also been noticed that the power low ([j]) describes only the general tendency of the function 
T{p — p c ), and one can observe some oscillations imposed on it. They can result from a ragged (fractal) measure of 
chaotic attractor colliding with its basin of attraction Q ; in the case of homoclinic crisis their period on a log-log plot 
of T(p — p c ) is I log IA2II, and their amplitude is large for small | A2 j , where A2 is the contracting eigenvalue of the crisis 
orbit. These oscillations arc indicated as a potential complication to verifying the scaling low (|l]) and determining 
the critical exponent 7. 

In this paper we investigate another kind of oscillations that, in general, can be more pronounced, and arises 
from an intertwined structure of the basins of attraction jio] [l2| . Their particular feature is the existence of sections 
where the average characteristic time T grows up when moving off the crisis point, opposite to the general trend 
([j]). Consequently, we call them anomalous oscillations. Oscillations of this kind have been observed in a ID map 
describing a diode resonator 0] and presumably in many other systems but, to our knowledge, they were not explicitly 
pointed out. We estimate the period and calculate the amplitude of these basin-induced oscillations using a simple 
model of an attractor colliding with its basin boundary. 



II. CRISIS IN A SPIN MODEL, ANOMALOUS OSCILLATIONS 

First let us briefly introduce a spin model in which the anomalous oscillations can be explicitly observed. After 
the papers |f^,|l4|], we consider a classical magnetic moment (spin) S, |S| = 5* in the field of uniaxial anisotropy (z— 
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easy/hard axis) with imposed transversal magnetic field B(t) along the x axis. The system can be described by the 
Hamiltonian H = — A(S Z ) 2 — B(t)S x , where A is the anisotropy constant. Such a model can describe a magnetic ion 
in a paramagnetic material or a single domain ferromagnetic sample. The motion of the spin is determined by the 
Landau-Lifschitz equation with damping term: 

— = SxB eff --Sx (SxB eff ), (2) 
at b 

where B e ff = —dH/dS is the effective magnetic field and A > is a damping parameter. 

Taking the driving field in the form of periodic delta pulses of amplitude B and period r: B(t) — B^2^ =1 5(t— nr), 
and using the fact that |S| is constant, the equation of motion (^) can be transformed into a superposition of two 2D 
maps [^^4|: Ta describing the time evolution between kicks and Tg describing the effect of the kick. The complete 
dynamics is yielded as a composition of the two maps: 

[g„<p]=T B [T A [S,Ml (3) 

where ip is the angle between the axis x and the projection of the spin on the xy plane. For different values of the 
parameters the system exhibits various types of dynamics including the periodic and chaotic ones |l5|]l4f| (see also 

As an example, consider the crisis that occurs at t c = 27r, A c = 0.1054942, A c = 1 and B c = 1 in which two 
symmetric chaotic attractors merge |L3|,|l4). The attractors correspond to two Ising states (spin "up" and "down") 
existing in the absence of the external field. We take the amplitude of the driving field B as an accessible system 
parameter. For B > B c random jumps between the two, previously separate attractors can be observed; 

Fig. [I] shows the dependence of the average time T between two subsequent switches on the distance B — B c from 
the crisis point in a log-log scale. We observe a linear trend according to ([[]) with the exponent 7 w 0.77 (dashed line) 
. However, a remarkable, roughly periodic oscillations around the trend line can also be seen. These oscillations 
include anomalous sections where T increases when moving away from the crisis point. 

The oscillations are due to an intertwined, fractal structure of the basins of attraction of the symmetric pre-crisis 
attractors. To see this, consider the situation depicted in Fig. || where a part of attractor is plotted together with 
pseudo-basin of attraction (gray spots) of the other, symmetric attractor. In fact, after crisis there is only one common 
attractor and almost all points of the phase space form its basin, but for B > B c we can define pseudo-basins as sets of 
initial conditions evolving to the "upper" (S z > 0) or the "lower" (S z < 0) attractor respectively after M iterations, 
where M <C T; in Fig. || we took M = 20. The structure of pseudo-basins is much similar to the real basins before 
crisis. In what follows, we shall refer to the post-crisis pseudo-basin of attraction of the other attractor as the basin 
of escape. The average transient time T is proportional to the inverse of the measure (i of the part of attractor 
overlapping with the basin of escape || . 

In Fig. U (a) the branch ^41 of the attractor crossed the band B3 of the basin of escape and is just before the next 
band B2. This corresponds to a local maximum in Fig. [j]. When the parameter B is increased Al enters the band 
B2, so \x grows and T decreases until Al touches the lower edge of B2 (Fig. |^ (b)). Then the total overlap /1 begins 
to decrease and we observe an anomalous section in Fig. ^. But simultaneously another branch A2 enters the basin 
of escape, so the anomalous section is not as pronounced as in the "clean" case when Al crosses B3. 

If the band S3 is magnified by a factor of 0.125 -1 a structure similar to that on Fig. || could be observed; the 
procedure could be repeated further on. This explains the fact that the mean period of oscillations in Fig. |l| is 
approximately | log 10 0.125| « 0.9. 



III. SIMPLE MODEL OF ANOMALOUS OSCILLATIONS 



Motivated by the observed self-similarity we introduce a simple model of an attractor creeping into an intertwined 
basin of escape in order to assess the amplitude of anomalous oscillations. We define the basin in the vicinity of a 
collision point as a self-similar set S of stripes of the width /3 l bE accumulating at the line y = (see Fig. ||): 

00 

» = Lkto y } '■ y > Ky< + ( 4 ) 

i=0 

P l bn is the width of the gaps in the basin of escape that may be called "return regions" . The parameters b, f3, bE and 
bji fulfil the condition bE + b^ + (3b = b. 
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As a model attractor we take a single parabola y — x 2 — r. The parameter r corresponds to the distance p — p c from 
the crisis point in ([!]). At r — the attractor collides with the basin 23 (crisis occurs). For r > there is a nonzero 
overlap of both sets. We assume the average transient time T(r) proportional to the inverse of the measure (i(r) of 
the overlapping region (cf. Fig. |] shows the function T(r) in log- log scale obtained from the described model. 
Due to self similarity of the model basin 23, we observe regular oscillations with the period log (3 superimposed on the 
linear trend described by 

Our aim is to calculate the amplitude A of the anomalous oscillations measured as the difference between a maximum 
and the preceding minimum, see Fig. g. Consider one fc-th period of the oscillation. By our definition 

A = logT(r 2 )-logT(r 1 ) = log4 !: 4, (5) 

where r 1 = (3 k b, a r 2 = /3 k b + f3 k - 1 b R (cf. Fig. g). 

For small r, the length of parabola in the half-plane y < is proportional to ^fr. Using this approximation the 
measure of the overlap fj, at points r% and r 2 can be written as 

A*(ri) ~ ^fWb~E, Mr 2 ) ~ y/^bR + ^bE - VP^r. (6) 
Here we neglected all the strips with i > k in ([|). From (g) and (g) we get 

A = (l/m^shm -1 y/RE//3, (7) 

where RE = 6r/6e and q is the base of the logarithm in the plot logT vs. logr (in this Letter we use q = 10). 

One can also calculate the slope of the linear trend, getting 7 = 5 because of a ID model-attractor. Analytical and 
numerical calculations considering other model-attractors show that the fractal structure of the basin of escape does 
not influence the slope 7. 



IV. NUMERICAL EXAMPLES: SPIN MAP AND HENON MAP 

The formula (Q) has been verified by comparing the calculated A with the amplitudes of anomalous oscillations 
measured on the log- log charts T(p — p c ), obtained from computer simulations of the spin map (g). We investigated 
the intermittent dynamics near four different crisis points of the same kind as the one described above, and we made 
respective plots T(B — B c ), similar to Fig. [l| The results are gathered in Table 1. The values of (3 and RE used in 
(0) have been measured employing a series of consecutive magnifications of the basin boundary at B = B c (similar to 
Fig. ||). Much the same values of /3 can be obtained as 10 - ^, where 9 is the average period of anomalous oscillations 
measured as in Fig. ^. 

Note, that the amplitudes of anomalous oscillations calculated from ((?]) are close to the maximal amplitudes 
observed in the simulations of the map (g) . This can be understood noting that the maximal amplitude corresponds 
to a "clear" situation when a branch of attractor which may be considered as one parabola (its fractal fine structure 
may be ignored) leaves a band of the basin while no other branch collides with any other basin band. It is just the 
case considered in our model crisis. Conversely, if one branch of attractor leaves a basin band and simultaneously 
another branch enters the same or some other band, both effects can cancel and no pronounced peak is observed. 
Thus, (0) gives the maximal amplitude of anomalous oscillations. 

We also studied a homoclinic boundary crisis that occurs in the Hcnon map, x n+ \ = p — x„ — Jy n ,y n +\ = x n 
at J c = 0.3 and p c = 2.124672450... |J||. If we make the log 10 -scale plot of T(p - p c ) (here T is the length of 
chaotic transient before the escape to infinity, averaged over a set of initial conditions) measuring T with appropriate 
accuracy, and marking subsequent points dense enough we can see tiny anomalous peaks with average period 0.97 
and average amplitude (defined as on Fig. [l]) 0.1. These oscillations are superimposed on the linear trend ([!]) and, 
here dominating, oscillations due to the ragged measure of the chaotic attractor. A few consecutive magnifications of 
the fractal basin boundary gives the model parameter values RE w 0.0167 and f3 ~ 0.107. The formula (g) then gives 
A = 0.168. This value is again close to the maximal amplitude observed (ss 0.16). Note, that the average period of 
anomalous oscillations is approximately | log 10 /3|. 



V. DISCUSSION AND SUMMARY 



The presented model used to obtain the maximal amplitude of anomalous oscillations is a simple one. In real 
basins of attraction every strip of 23 may have its own fine structure and chaotic attractor is locally an infinite fractal 
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set of parabolas. Nevertheless, the model comprises the main properties of both sets responsible for the anomalous 
oscillations and gives a proper estimate of their maximal amplitude. 

Anomalous oscillations seem to be quite common in nonlinear dynamical systems, as do the fractal, intertwined 
basin boundaries that give rise to them. They are, however, easy to confuse with statistical dispersion of the mean 
times T, especially when one has a small number of points on T(j> — p c ) plots. The presence of anomalous oscillations 
can be a complication to verify the scaling low (Q), but if one is able to detect and measure them, they can give 
additional information on the structure of basins of attraction near crises. Note, that both A and (3 can be, in 
principle, assessed from the plot (like in Fig. [I]) and Eqn. (0) can be used to derive RE. 

A model of chaotic attractor similar to that of the basin of attraction lets us calculate the amplitude of "normal" 
oscillations caused by the fractal structure of the attractor [|l6| . 

Concluding, we have observed anomalous oscillations of average characteristic crisis-induced transient times, and 
explained them in terms of a simple, geometric model, which gives the period of the oscillations and enables to 
calculate with a fair accuracy their maximal amplitude. The model is universal, and the results are applicable to a 
large class of 2D systems undergoing different types of crises. It could be also customized for some special non generic 
cases. 
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FIG. 1. The average time of residence on one of the pre-crisis attractors versus the distance from the crisis point 
A c = 1,T C = 2ty, A c = 0.1054942..., B c — 1 for the spin map (p|). Standard deviation is of order or below the size of the 
plotted points. Two parameter values which correspond to Fig. Hare marked. 
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FIG. 2. Attractor of the spin map (|3|) creeping into the pseudo-basin of attraction of the other, symmetric attractor at 
A = l,r = 2tt, A = 0, 1054942... and (a) - B = 1,0000032, (b) - B = 1,000015. Am (m = 1,2,...) denotes branches of the 
attractor, and Bm bands of the basin of attraction. 



FIG. 3. Model of an intertwined basin boundary. 

FIG. 4. Function T(r) obtained when a model one-parabola attractor enters the model pseudo-basin of attraction 93 after 
crisis. 





Crisis parameters 




Model parameters 


Amplitude of oscillations 




B 


T 


A 


P 


RE 








1 


2tt 


0.10549 


0.125 


0.5 


0.502 


0.63 


0.627 


0.92 


2tt 


0.08964 


0.08 


0.95 


0.7 


0.8 


0.847 


0.92 


5.8866 


0.08 


0.16 


0.15 


0.325 


0.37 


0.373 


1.2 


2tt 


0.1437 


0.124 


0.077 


0.26 


0.31 


0.314 



TABLE I. Amplitudes of anomalous oscillations of average characteristic times T in the vicinity of four different crisis points 
in the spin map (pi) at A — 1. The average (A a „) and maximal (A max ) amplitudes are compared to the values Ath obtained 
from (T7J) with respective (3 and RE. 
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